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Solar wind test of the de Broglie-Proca's massive photon with Cluster 

multi-spacecraft data 
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We use spacecraft data in the solar wind at 1 AU to provide an estimate of the upper limit of 
the de Broglie-Proca's massive photon, by looking for deviations from the classical Ampere's law. 
We take advantage of the Cluster spacecraft which allow the direct computation of V x B from 
simultaneous four-point measurements of the magnetic field. We estimate the upper bound for the 



mass of the de Broglie-Proca's photon m 7 to be 2 ■ 10 
based on ad hoc models in the solar wind. 



kg, in agreement with previous findings 



PACS numbers: 03.50.De, 12.20.Fv, 14.70.Bh, 95.55.-n, 96.50.Bh, 96.50.Ci 



Our understanding of the universe is largely based on 
electromagnetic observations. As photons are the main 
messengers, fundamental physics has a concern in testing 
the foundations of clcctromagnctism. Non-Maxwellian 
theories can be grouped into two main classes: non-linear 
and massive photon theories. The former constitute a 
branch of classical and quantum field theories, and from 
their conception by Born and Infeld [l[, Heisenberg and 
Euler 0] , non-linear theories have been at the heart of the 
description of a plethora of phenomena. The latter class 
assume or predict the photon being massive, whereas for 
the Standard Model the only massless particles are the 
photon and the graviton, i.e. the bosons of the eldest 
discovered interactions. 

The concept of a massive photon has been vigorously 
pursued by Louis de Broglic from 1922 Q throughout 
his life. He defines the value of the mass to be lower 
than 10 -53 kg A comprehensive work of 1940 [|[ 
contains the modified Maxwell's equations and the re- 
lated Lagrangian. Instead, the original aim of Alcxandru 
Proca [37| . was the description of electrons and positrons. 
Despite Proca's several assertions on the photons be- 
ing massless, his Lagrangian @ and formalism Q ap- 
ply to a massive real or complex vector field. Theories 
and conjectures centered on massive photons have been 
later proposed by several authors. A non-exhaustive list 
of gauge invariant formalisms includes that of Podolsky, 
who discuss higher derivatives of the gauge field [8|; of 
Stueckelberg, who insures gauge invariance by adding a 
scalar field @ ; of Chern and Simons [l(| , whose elec- 



trodynamics couples the field and the potential in the 
Lagrangian through the Levi-Civita's tensor. Recent re- 
views can be found in [llHIl. The concept of massive 
photons impacts on many fields of physics and astro- 
physics, e.g., charge conservation and quantization, mag- 
netic monopolcs, superconductors, charged black holes 
and the cosmic microwave background. The relation be- 
tween the massive photon and the Higgs' boson has been 
also addressed [lU, [l5| . The first experimental limits to 
photon mass were provided by Schrodinger, who pointed 
out the finiteness of the range of the modified electro- 
magnetism (l6l 

How much the theory of relativity would be affected by 
the massive photon assumption, is not straightforward to 
assess: partly due to the variety of the theories above, 
and partly to the removal of our ordinary landmarks and 
the rising of interwoven implications. Indeed, a phe- 
nomenologically oriented survey displaying the testable 
differences among different theories is missing, and exper- 
imentalists have mostly conveyed their efforts in checking 
the simplest and eldest massive photon model, i.e. the 
photon by de Broglie-Proca, henceforth dBP. The upper 
limits to the photon mass found in such experiments can- 
not be generalized to other massive photon theories. The 
massive electromagnetism dBP equations for the electric 
E and magnetic B fields were written first by de Broglic 
@. In SI units, they are 



V -E = i- -M 2 4> , 



(la) 



'Email: alcssandro.rctino@lpp.polytcchniquc.fr 
Web-page: http:/ /www.lpp.fr/?Alessandro-Rctino 
i Email: spallicci@cnrs-orleans.fr 
Web-page: http: / /lpc2e. cnrs-orleans.fr/~spallicci/ 



V • B = , 



(lb) 



(lc) 



2 



VxB = fi j + /i eo 



BE 
~dt 



MA 



(Id) 



where no is the free space permeability, p the charge den- 
sity, j the current vector, 4> and A are the scalar and vec- 
tor potentials, M. = m^c/h = 1/A, h being the reduced 
Planck's constant, c the speed of light, A the reduced 
Compton's wavelength, and m 7 the photon mass. Equa- 
tions (|laHld|) are Lorentz but not gauge invariant, due 
to the explicit presence of the potentials. This implies 
that the gauge must be fixed when performing a mea- 
surement. Two laws are modified with respect to the 
original Maxwell's formulation: the curl of the magnetic 
field (Ampere, 1826 [3; Maxwell, 1861 [H), and the 
divergence of the electric field (Gauss, discovered in 1835 
and published in 1867.) The dates of these 19 th century 
achievements are in striking contrast with the current 
complex and multi-parameterized cosmology. 

Experiments have constrained the photon mass to very 
low limits. A commonly accepted limit from a laboratory 
experiment (Coulomb's law) is 2 x 1CP 50 kg [20] . One 
prediction of the dB-P theory affirms that lower energy 
photons travel at lower speed than those at higher en- 
ergies. This prediction behaves like the dispersion from 
plasma. In pulsar timing, different arrival times of the in- 
coming photons are routinely measured, but lacking any 
other independent measurement on the electron density, 
the differences are solely attributed to plasma dispersion. 
The dispersion-based limit is 3 • 10 -49 kg [2l|. A num- 
ber of different experiments and analysis have been per- 
formed on space environments and obtained lower lim- 
its. Davis et al. [22[ use Pioneer-10 data of the Jupiter 
magnetic field to set a limit of 8 • 10~ 52 kg. Tighter 
estimates emerge by the analysis of the interplanetary 
magnetic field in the solar system, although Adclbcrgcr 
et al. [lH question such limits if the photon were to ac- 
quire its mass through the Higgs' mechanism. Ryutov 
found first m 7 < 10 -52 kg in the solar wind at 1 AU 
[H, and later m 7 < 1.5 x 1CT 54 kg at 40 AU [H. The 
latter is the limit currently accepted by the Particle Data 
Group (PDG) (25)]. Even lower limits have been claimed 
from modeling the galactic magnetic field and are 10~ 62 
kg [TH [2^1 . It should be noted that the lowest theoretical 
limit on the measurement of any mass is prescribed by 
the Heisenberg's uncertainty principle m > h/ Ate 2 , and 
gives 2.7 x 10~ 70 kg, where At is the currently supposed 
age of the Universe (1,37 • 10 10 years). Incidentally, the 
same principle implies that measurements of masses in 
the order of 10 -54 kg should be performed in time scales 
of about twenty minutes. 

An important warning must be flagged. The exami- 
nation of the literature on the large-scale " astrophysical 
limits" inspires a critical attitude and prompt to question 
whether these limits are nothing more than the outcome 
of idealized models. This view is confirmed by Goldhabcr 
and Nieto I38| . Estimates from planetary [22|, [27| and so- 
lar wind |23l. |24T| magnetic fields are likely to be more 



reliable, as they are based on in situ measurements and 
the estimated upper limit of p hoton mass m 7 is closer 
to laboratory experiments [20|. But even the estimates 
from in situ measurements are based on several assump- 
tions. For the case of the solar wind, a close scrutiny 
of the last accepted limits [23|, [24[ reveals that: (i) the 
magnetic field is assumed exactly always and everywhere 
a Parker's spiral; (ii) the accuracy of particle data mea- 
surements (from e.g. Pioneer or Voyager) has not been 
discussed; (iii) there is no error analysis. It is then impor- 
tant to verify the robustness of such solar wind estimates 
by using a more detailed experimental analysis. 

Dealing with such a small mass should induce to ex- 
treme caution. A small mass needs a very precise exper- 
iment or, alternatively, a very large apparatus, since a 
small mass is associated to a very large (reduced) Comp- 
ton wavelength A. In this letter, we focus on the second 
possibility through the largest-scale magnetic field acces- 
sible to in situ spacecraft measurements, i.e. the inter- 
planetary magnetic field carried by the solar wind. For 
this purpose, we use the only currently available multi- 
point spacecraft to test the dBP modified Ampere's law. 
Cluster [28j is an European Space Agency mission com- 
posed by four identical spacecraft flying in tctrahedral 
configuration. Cluster has allowed for the first time the 
direct estimate of three-dimensional quantities from four- 
point in situ measurements of particle and electromag- 
netic fields [28| . The novelty of this approach is that we 
can directly estimate V x B from magnetic field mea- 
surements; this was not possible with earlier spacecraft 
measurements. For the steady components of the mag- 
netic field, i.e. low frequencies, the displacement current 
term can be dropped, and we thus refer to the Ampere's 
law. Its dBP modified version reads 



V x B = n Q j - M 2 A 



(2) 



The only experiment able to properly estimate from 
Eq. [2] the mass m 7 requires the independent, and within 
experimental errors, measurement of the difference be- 
tween the currents js = V X B//io and j = jp = 
ne(vi —v e ), where n is the number density, e the electron 
charge, and Vi, v e are the velocity of the ions and elec- 
trons, respectively. To our knowledge, this measurement 
has never been performed up to now. The dBP photon 
mass is given by (SI units) 




ne(vi - v e ) 



V x B 



/'o 



JP - JB 



(3) 

where h = 1.05 • 10~ 34 J s, c = 2.99 • 10 8 m s" 1 , = 
1.26- 10 -6 N A -2 , A H is the vector potential of the inter- 
planetary magnetic field, and k = ft/i^c -1 \ Ah \^ 1/ " 2 . 

The event selection has been performed according to 
the following criteria: (i) a calm solar wind event, i.e. 
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FIG. 1: Panel (a). Three components of the magnetic field for 
Cluster 1 in the GSE (Geocentric Solar Ecliptic) coordinate 
system. The magnetic field is measured at 22 samples/s and 
has been re-sampled at the time resolution of the spacecraft 
position measurements (1 minute). The GSE coordinate sys- 
tem is defined as having the X direction pointing from the 
Earth to the Sun, the Z direction orthogonal to the ecliptic 
plane, and the Y direction completing the right-handed sys- 
tem. Panel (b). The B x component (solid line), and the same 
component low-pass filtered below 0.01 Hz (dashed line). 



disconnected from the terrestrial bow shock and as far 
as possible from the terrestrial magnetic field; (ii) the 
location of the four spacecraft closest to the equatorial 
plane; (iii) the widest inter-spacecraft separation, namely 
10 4 km, for the largest differences in the magnetic field 
among spacecraft; (iv) the configuration best approach- 
ing the tetrahedron; (v) the availability of the ion and 
electron moments. Figure [T] shows the magnetic field 
components measured by Cluster for one such event dur- 
ing a time interval of two hours, in the in GSE (Geocen- 
tric Solar Ecliptic) coordinate system. It is to be inferred 
from Fig. [TJ that the magnetic field has B x < 0, B y > 
and B x ,B y ^> B z , as expected for the Parker's spiral 
close to the ecliptic plane, where Cluster passes at the se- 
lected event. It is emphasized, however, that our analysis 
does not rely on the Parker's model, since the magnetic 
field is measured in situ. The choice of one event having 
Parker's spiral orientation is to have conditions as similar 
as possible to those present in [23|, [24[. For this event, 
A H ~< | B 1 > xL = 2.65 • 1CT 9 x 1.5 • 10 11 = 4 • 10 2 
T m, where < |£?o| > is the value of the magnetic field 
measured by Cluster, Fig. [T] and L, distance of Cluster 
from the Sun, is 1 AU in meters. 

In this paper we estimate V x B by using the curlome- 
ter method [23, |3(| on the magnetic field from the Cluster 
flux-gate magnetometer data [3l| . This method allows (i) 
to compute the average Vxfi over the spacecraft tetra- 
hedron with no assumptions on the field analytical form 
(only assuming linear gradients); (ii) to assess the error 
on the average of the current V x Bj fj,Q 
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FIG. 2: Panel (a). The current density computed with the 
curlometer method. The dashed horizontal line corresponds 
to its average value. Panel (b). Estimate of the deviations 
from linearity for the curlometer method, where the lower 
values are associated to higher linearity. Panel (c). Rela- 
tive error for the current of the current density. Panel (d). 
Elongation of the tetrahedron < E < 1. An elongation of 
zero, corresponds to a perfect tetrahedron. Panel (e). Pla- 
narity of the tetrahedron < P < 1. Planarity equals to zero 
corresponds to a perfect tetrahedron. 



with AB being the error on the average of the magnetic 
field, and AR the error on the average of the spacecraft 
separation R; (iii) to estimate the deviation from the lin- 
earity assumption through the quantity |V • -B|/|V x B\. 
Figure [2] shows the result of the curlometer for the 
event of Fig. [TJ after having low-pass filtered the mag- 
netic field with fip = 0.01 Hz to remove the magneto- 
hydrodynamic (MHD) waves, as we are interested to the 
steady components of the magnetic field only. 

Our appraisal of j B = \VxB\/fi is 1.37- 10- u ± 5.27- 
10 -13 A m -2 . The relative error Ajs/ < 3b > is about 
5%, obtained by assuming AB = 0.1 nT (accuracy of the 
magnetometer), < R >= 10 4 km (spacecraft separation), 
AR = 1% x < R >= 100 km. The accuracy of this esti- 
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mate is satisfactory, as indicated by the small relative er- 
ror, and by the reasonable value of | V-£?|/| V x B\ ~ 50%, 
in agreement with earlier accuracy estimates |30j . The 
values elongation ~ 0.2 planarity ~ 0.5 of the tetrahe- 
dron are within the range for accurate estimates (30j |. 

Conversely, a meaningful assessment of the current 
density from the ion and electron moments is not pos- 
sible due to inherent instrument limitations. Ion electro- 
static analyzers, such as that on-board Cluster [23], can 
saturate in the solar wind due to the presence of high 
fluxes, and thereby resulting often in incorrect ion mo- 
ments. Further, the computation of the moments from 
the electron electrostatic analyzers on-board Cluster is 
affected by photoelectrons and spacecraft charging issues 
[33j . The assessed value of the current from the parti- 
cle detectors for this event (as well as for other events) 
is much larger than those from the curlometer (this is 
mostly due to the differences in the velocities, while the 
estimate of particle density is reasonable). These large 
values are not consistent with the values expected from 
the difference | jp — jp |, based on earlier findings of 
the upper limit of to 7 , in the range 10 -54 - 10~ 50 kg. 
Besides, measurements in other regions of space, e.g. 
Earth's magnetosphere, typically show that this differ- 
ence is small [13| • We have therefore disregarded particle 
instrument data to avoid unjustified claims on experi- 
mental limits, or even more, on the discovery of a massive 
photon. Nevertheless, we can still put a conservative up- 
per limit on m 7 . We write (for A indicating the absolute 
error) 

m 7 < k [j P + Aj P - {j B - Aj B )]* . (5) 
We assume that jp — jp is very small, and thus 

m i < k [is (ep + ep)] 2 , (6) 

where ep.p arc the relative errors. We further assume 
that such errors are comparable. A worst case estimate 
based on particle data leads to ep ~ 1, and thereby im- 
plying a factor 3 on the upper limit of m 1 . Relying on 
these two assumptions, the estimate comes from the error 
of the curlometer 

m 1 < k(2e B j B )i = 2.83 ■ 10" 47 (Aj B )* kg, (7) 

that is m 7 < 2 • 10~ 53 kg, for Ajp = 5.27- 10~ 13 A m" 2 . 

This value is one order of magnitude smaller than the 
earlier estimate in the solar wind at 1 AU [23l | . but one 
order of magnitude larger than the estimation done at 40 
AU [13], that is the currently accepted limit by PDG [25j |. 
Our value of m 7 from a detailed experimental analysis is 
thus in agreement with such estimates coming from ad 
hoc models. On the other hand, our estimate is three 
orders of magnitude smaller than laboratory limit [2Cj |. 

Improvements of this estimate could come from a bet- 
ter accuracy of spacecraft separation and, most impor- 
tantly, of magnetic field measurements, as well as from 



more accurate particle measurements. We have consid- 
ered a conservative error for the spacecraft separation of 
1%, corresponding to about 100 km. Nevertheless, typ- 
ical values from the literature are smaller, i.e. 5 — 10 
km [30| . thereby reducing the contribution of the space- 
craft position error to Ajp by a factor of 10 — 20. How- 
ever, since the largest contribution to Ajp comes from 
the error on the magnetic field measurements, the im- 
provement in the accuracy of separation would not be 
crucial. A better accuracy of the measurements from the 
fluxgate magnetometer is more consequential. For our es- 
timate, we have conservatively adopted AB = 0.1 nT as 
typical accuracy for the fluxgate magnetometer on-board 
Cluster. An accuracy of 0.01 nT would result in an im- 
provement factor of about 3 on the upper limit of m 7 . 
Such accuracy is already included in the specifications of 
fluxgate magnetometers, such as those on-board Cluster 
[3l| . or Cassini [35)] spacecraft. Further advances may be 
achieved in the future by more accurate magnetometers. 

Undoubtedly, the most stringent requirement comes 
from the particle detectors. For the event studied here, 
the difference between ion and electron velocities is 



\vi-v e \ = J —~ 6.2 ■ 10 11 jp m s-\ (8) 

ne 

where n ~ 10 7 m 3 (ion density measured for this event) 
and e = 1.6 • 10~ 19 C. By assuming the current density 
to be similar to that given by the curlometer, that is 
jp = jp = 1.37 • 10~ n A m -2 , the difference of ion and 
electron velocities should be less than 10 m s — 1 . This 
velocity difference cannot be measured neither by the 
electrostatic analyzers on-board Cluster, nor by more 
recent ones. The typical energy resolution 5E/E of 
such detectors is several percents [H, [33|j. By assum- 
ing 5E/E ~ 10%, and a typical value of particle energy 
in the solar wind E ~ 10 eV, the resolution for the ve- 
locity is Sv e <~ 95 km s -1 for the electrons, and Svi ~ 2 
km s _1 for the ions. It is straightforward to see that such 
detectors cannot resolve a velocity difference of 10 m s _1 . 

Few other considerations are worth mentioning. The 
Planck and the permeability (or magnetic) constants are 
depending on the definition of the speed of light which in 
itself depends on the photon mass. Therefore, there is an 
underlying epistemological contradiction in the equations 
above. Given the size of the errors in the currents and 
in the vector potential, the error contributions from the 
constants were neglected. 

Most non-Maxwellian theories foresee a deviation from 
the Ampere's law. An analysis of these theories versus 
Cluster data is planned [39[ . These investigations empha- 
size the importance of a dedicated space measurement. 

We have reported a novel test for estimating the up- 
per limit of the dBP massive photon by using Cluster 
multi-spacecraft measurements in the solar wind. We 
have found to 7 < 2 ■ 10 -53 kg, in agreement with pre- 
vious solar wind estimates, but our test being based on 
fewer assumptions than previous results. First, we have 
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directly assessed V xB/fio from four-point magnetic field 
measurements; this was not possible before. Second, our 
test does not assume that the interplanetary magnetic 
field is a Parker's spiral. Third, is the only measurement 
performed so far in the solar wind that takes into account 
in detail the effect of experimental errors. Finally, though 
the de Broglic's prediction on the upper limit of the pho- 
ton mass is confirmed [§| , the domain between the labora- 
tory findings (m 7 < 10 -50 kg) and the solar wind results 
(to 7 < 10~ 54 kg) is still subjected to assumptions and 
conjectures, though fewer now, and not to clear-cutting 



outcomes from experiments. Our experiment is limited 
by the resolution of the velocity differences of ions and 
electrons. 
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